Dynamic polarizabilities for the low lying states of Ca"*" 
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The dynamic polarizabilities of the 4s, 3d and 4p states of Ca"'', are calculated using a relativistic 
structure model. The wavelengths at which the Stark shifts between different pairs of transitions 
are zero are computed. Experimental determination of the magic wavelengths can be used to 
estimate the ratio of the f^dj^Apj, and /i^^ .^-^ip , oscillator strengths. This could prove valuable 
in developing better atomic structure models and in particular lead to improved values of the 
polarizabilities needed in the evaluation of the blackbody radiation shift of the Ca^ ion. 

PACS numbers: 31.15.ac, 31.15.ap, 34.20.Cf 



I. INTRODUCTION 



The dynamic polarizability of an atom or ion gives a 
measure of the energy shift of the atom or ion when im- 
mersed in an electromagnetic field [ll-l3l|. For any given 
state, one can write 



AS = --arfHF2 



(1) 



where adi^o) is the polarizability of the quantum state 
at frequency w, and F is a measure of the strength of 
the AC electromagnetic field. The value of the dynamic 
polarizability in the a; — > limit is the static dipole po- 
larizability. 

The magic wavelength for a transition is the wave- 
length for which the AC Stark shift of the transition 
energy is zero [^-u|. The identification of magic wave- 
lengths and their use in making optical lattices has re- 
sulted in the development of optical lattice clocks which 
have the potential to exceed the performance characteris- 
tics of the existing standard for time, namely the cesium 
microwave clock |8l-[l2|. 

However, the experimental determination of magic 
wavelengths also provides valuable information to con- 
strain the atomic structure models that are used to es- 
timate the impact of Stark shifts on the performance on 
atomic and ion clocks [Ij, [lj| . A parameter related to 
the magic wavelength is the tune-out wavelength. The 
tune-out wavelengths for an atomic state are the wave- 
lengths at which the polarizability for that state goes to 
zero [15-17]. It should be noted that most atomic states 
have a number of tune-out wavelengths just like most 
atomic transitions have a variety of magic wavelengths. 

The advantage of magic and tune-out wavelength mea- 
surements are that they are effectively null experiments. 
They measure the frequencies at which polarizability re- 
lated quantities are equal to zero. Therefore they do not 
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rely on a precise determination of the strength of a static 
electric field or the intensity of a laser field. This makes 
it possible to determine the magic wavelengths to a high 
degree of precision d, [H-llll. 

There have been a number of theoretical studies of 
the properties of the low-lying Ca"*" ion |23l- [29[ by 3 dif- 
ferent research groups. One of these groups 24| used a 
non-relativistic approach while the other two groups used 
explicitly relativistic formulations |26l - l29l |. One of the 
singular features about the relativistic calculations are 
significant differences between predictions of the prop- 
erties of spin-orbit doublets. The relativistic all-order 
many body perturbation theory method predicts rela- 
tively small non-geometric differences between the line 
strengths of the idj and Apj spin-orbit doublets [23, |23] . 
The relativistic coupled cluster approach typically gives 
much larger differences [26[ . One of the secondary aims 
of the present work is to shed light on these differences. 

The present manuscript reports calculations of the dy- 
namic polarizabilities of the five lowest states of Ca+. 
The Hamiltonian used is a fully relativistic version of 
a semi-empirical fixed core potential that has been suc- 
cessfully applied to the description of many one and 
two electron atoms |30l - [33| . While there are many dif- 
ferences in the technical detail, the underlying philoso- 
phy and the effective Hamiltonian for the valence elec- 
tron are essentially the same once the relativistic mod- 
ifications are taken into account. Magic wavelengths 
for the 4s — ?> 8^3/2,5/2 and 4s — ;> '^Pi/2,3/2 transitions 
are given. The dynamic polarizability of the ground 
Ca"'"(4s) state is dominated by the 4s -^ Apj transitions 
and its accuracy is largely dependent on the accuracy 
of the transition matrix elements connecting the 4s and 
4pj states. The description of transitions involving the 
Ca+(3(i) state is complicated by the effect that the 3d 
electrons have on the core electrons. The 3d orbitals have 
the smallest (r) expectation values of any of the valence 
electrons and this does distort the wavefunctions for the 
outermost core electrons 34, 35j. One consequence of 
this are greater uncertainties in the calculation of transi- 
tion matrix elements involving the 3dj states [23, [23, |33 



All results reported in this paper are given in atomic 
units with the exception of the lifetimes which are given 
in seconds. The value adopted for the speed of light is 
c = 137.035999074(44) a.u. 



II. FORMULATION AND ENERGIES 

A. Solution of the Dirac-Fock equation for closed 
shell atomic system 

The calculation methodology is as follows. The first 
step involves a Dirac-Fock (DF) calculation of the Ca^+ 
ground state. The DF calculation begins with the equa- 
tion 

' N ^ ^ 

^i/D. + V— UM-i^V'M, (2) 






i<j 



«J 



where Hdi is the single-electron Dirac Hamiltonian 

Hj,.^ = ca, ■ p, + c2(/3, - 1) + Vin). (3) 

where c is the speed of light, p is the momentum operator, 
and a and /? are the Dirac matrices pBg . 

The orbitals of the DF wave function, i!{r), can be 
written as 



V'(r) 






(4) 



where gnuir) and /««(?") are the large and small compo- 
nents, ^Km{f) and n_Km(J^) correspond to the angular 
components. The radial Dirac equation for an orbital 
can be expressed schematically as 



'V{r) + VMr) -<fr--r) 

^ c(f + ^) ^2c^ + V{r) + VMr) 



' gnAr) ' 



where Vdf is called the Dirac-Fock potential, and V{r) 
is the interaction potential between the electron and the 
nucleus. A Fermi nuclear distribution approximation is 
usually adopted for many-electron atomic system. 

The single particle orbitals are written as linear combi- 
nations of analytic basis functions and so the method of 
Roothaan [33, [33] is used to recast the DF equations into 
a set of matrix equations. The functions chosen are B- 
splines with Notre-Dame boundary conditions [39|. The 
large and small components are expanded in terms of 
a B-spline basis of k order defined on the finite cavity 

U , -Tlniax 1 



N 



^,{r)^J2^t''B.Mr) 



funir) 



N 
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CrB^Ar)■ 



(6) 



(7) 



TABLE I: Theoretical and experimental energy levels (in 
Hartree) for some of the low-lying states of Ca"*". The en- 
ergies are given relative to the energy of the Ca^""" core. The 
experimental data were taken from the NIST tabulation [43|. 



Level 


DF 


DFCP 


Experiment [43] 


4Si/2 


-0.4166315 


-0.4362777 


-0.4362776 


3^3/2 


-0.3308695 


-0.3740834 


-0.3740827 


3^5/2 


-0.3307597 


-0.3738074 


-0.3738062 


4pi/2 


-0.3099986 


-0.3214966 


-0.3214966 


4p3/2 


-0.3090889 


-0.3204818 


-0.3204810 


5Si/2 


-0.1933158 


-0.1983486 


-0.1985876 


4^3/2 


-0.1687383 


-0.1751536 


-0.1772989 


4d5/2 


-0.1686641 


-0.1750622 


-0.1772114 


5pi/2 


-0.1567656 


-0.1603178 


-0.1604688 


5p3/2 


-0.1564329 


-0.1600612 


-0.1601123 



The finite cavity is set as a knots sequence, {ti}, satis- 
fying an exponential distribution [40, |4l| . The specifics 
of the grid were that i?inax = 60 Oq and 50 B-splines of 
order k = 7 were used to represent the single particle 
states. Using the Galerkin method and MIT-bag-model 
boundary conditions [33| , the DF equations were solve by 
iteration until self-consistency was achieved. The single- 
electron orbital (Koopmans) energies of the closed shell 
Ca^"*" ion agreed with those computed with the GRASP92 
program [44I to better than 10~^ a.u. 



B. Polarization potential 

The effective potential of the valence electron with the 
core is then written 



Vc, 



Vkir(r) + yoxc(r) + Vpoi(r) 



(8) 



?he /direct and exchange interactions of the valence 
cKctVoii I ■«i(3;)r the DF core were calculated exactly. 
■Tlie /-dependent polarization potential, Vpoi, was semi- 
empirical in nature with the functional form 






ogl{r) 



2r4 



mm- 



(9) 



The coefficient, acorc is the static dipole polarizability 
of the core and gjAr) = 1 — exp(— r^/p? ) is a cutoff 
function designed to make the polarization potential fi- 
nite at the origin. The static dipole polarizability core 
was set to acorc = 3.26 a.u. [23. The cutoff parame- 
ters, pej were tuned to reproduce the binding energies 
of the ns ground state and the npj^ ndj excited states. 
Values of the cutoff parameters are po,i/2 — 1.7419 ao, 
PiA/2 = 1-6389 ao, Pi^i/2 = 1.6354 ao, P2,3/2 = 1-8472 ao, 
and P2,3/2 — 1-8489 ao. The cutoff parameters for ^ > 3 
were set to a common values of 1.897 qq. Table |T] gives 
the calculated B-spline and experimental energies com- 
ing from [43| . The calculations with the core-polarization 
potential are identified as the Dirac-Fock plus core polar- 
ization (DFCP) model. Differences between DFCP and 



TABLE II: Comparison of the electric dipole (El), electric 
quadrupole (E2) reduced matrix elements of several interested 
states of the Ca""" ion. 



Transition 



DFCP 



MBPT-SD 



RCC 



4Si/2 
4Si/2 
4Sl/2 

4si/2 ■ 



3d; 
3d 



'3/2 
'3/2 ■ 



3d5/2 
3d3/2 



3d; 
3d 
3d 



'3/2 ■ 
'5/2 ■ 
'3/2 ■ 



3d5/2 
3d5/2 
4pi/2 

4pi/2 

4p3/2 
4p3/2 



■ 4pi/2 

■ 4p3/2 
5pi/2 
5p3/2 

■ 4pi/2 
' 4p3/2 

- 4p3/2 

- 5pi/2 

■ 5P3/2 

■ 5p3/2 

■ 4/5/2 
' 4/5/2 
' 4/7/2 

- 5Si/2 

3d. 



4d; 
4di 



'3/2 
'3/2 
5/2 



4Sl/2 - 3d3/2 
4si/2 - 3d5/2 



2.879 
4.073 
0.089 
0.109 
2.500 
1.116 
3.356 
0.091 
0.044 
0.131 
1.964 
0.526 
2.354 
2.081 
4.205 
1.894 
5.675 
( 
8.120 



Dipole 

2.898(13) [29] 
4.099(18) [29] 



2.464(16) [29] 
1.100(6) [29] 
3.306(18) [29] 



1.927(52) [29] 

0.516(6) [29| 

2.309(29) [29] 

2.073(11) [29] 

4.28(3) [29] 

1.93(1) [29] 

5.78(3) [29] 

luadrupole 

7.939(37) [23] 



2.88(1) [26] 
4.03(1) [26] 



2.40(2) [26] 
1.09(1) [26| 
3.22(4) [26] 



9.964 9.740(47) [23] 



7.973 [44] 
8.12(5) [45] 

9.979 [44] 
9.97(6) [45] 



experimental energies mostly occur in the fourth digit 
after the decimal point. 

One of the interesting aspects of Table H] concerns the 
spin-orbit splitting of the ipj and 5pj states. The polar- 
ization potential parameters Pi,i/2 and pi,3/2 were tuned 
to give the correct spin-orbit splitting of the Apj states. 
Making this choice resulted in the spin-orbit splittings 
for the 5pj states also being very close to experiment. 



III. TRANSITION MATRIX ELEMENTS AND 
ASSOCIATED QUANTITIES 

A. Reduced Matrix Elements 

The dipole matrix elements were computed with a 
modified transition operator [30, H, SJ , e.g. 



rC^ ^rC^ - (1 



exp(-r7p'')) 5 



(10) 



The cutoff parameter, p used in Eq. ([T0| was set to p = 
ipta,ja + Pibdh)/"^ where a, b refer to the initial and final 
states of the transition. 

The static quadrupole polarizability of the Ca^+ core 
is needed for the calculation of the lifetimes of the 3dj 
states. It was set ag_corc — 6.936 a.u. |48j . 

There have been a number of previous calculations of 
reduced matrix elements and polarizabilities for the low- 
lying states of Ca"*" . The semi-empirical configuration in- 
teraction plus core polarization (CICP) can be regarded 



TABLE III: Comparison of the line strengths ratios for tran- 
sitions involving various spin-orbit doublets. The notation 
4si/2 — 4p3/2:i/2 means the line strength ratio defined by di- 
viding 4p3/2 line strength by the 4pi/2 line strength. 



Transition 



DFCP 



MBPT-SD 



RCC 



4Si/2 - 4p3/2:l/2 
4Sl/2 - 5p3/2:l/2 
3d3/2 — 4p3/2:l/2 



4p3/2 
4Si/2 

4P3/2 
4p3/2 



3di 
3d: 



'5/2:3/2 
'5/2:3/2 



4di 
4di 



'5/2:3/2 
'5/2 



2.0014 
1.4990 
5.0182 
9.043 
1.5057 

5.0181 
9.021 



2.001 [29] 1.958(17) [26] 



5.02 [29] 

9.04 [29] 

1.5052 [29] 



4.85(12) [26] 

8.73(27) [2^ 

1.5665 [44j 

1.5075 [45] 



as a non-relativistic predecessor of the present calcula- 
tion [2J, |30|- Another method used is the relativistic 
all-order single-double method where all single and dou- 
ble excitations of the Dirac-Fock (DF) wave function are 
included to all orders of many-body perturbation theory 
(MBPT-SD) [Uli^,!!!]. There have also been calcula- 
tions using the relativistic coupled cluster (RCC) method 
[H. The RCC and MBPT-SD approaches have many 
common features [50l - [53 |. Atomic parameters computed 
using the RCC approach have on a number of occasions 
had sig nificant differences with independent calculations 

[Uliilii]. 

The reduced matrix elements between the various low 
lying states are the dominant contributor to the polar- 
izabilities of the 4s, 3d and 4p levels. These are given 
in Table |TT] and compared with the results from other 
recent calculations. The ratio of line strengths for spin- 
orbit doublets is also interesting to tabulate since they 
can reveal the extent to which dynamical effects (as op- 
posed to geometric effects caused by the different angular 
momenta) are affecting the matrix elements. Some line 
strength ratios are given in Table IIIIl 

The variation between the DFCP, MBPT-SD and RCC 
matrix elements listed in Table Hill does not exceed 5%. 
The DFCP matrix elements are usually closer to the 
MBPT-SD calculations than the RCC matrix elements. 
A better indication of the differences between the DFCP, 
MBPT-SD and RCC calculations is gained by examina- 
tion of the line strength ratios listed in Table IIIIl The 
DFCP line strength ratios are within 1% of the values 
that would be expected simply due to the angular mo- 
mentum factors alone. The ratios are in very good agree- 
ment with the MBPT-SD ratios. It should be noted, 
that the line strength ratios for the resonant transition 
of potassium have been measured to be very close to 2.0 
[13 and DFCP and MBPT-SD calculations also predict 
line strength ratios very close to 2.0 |17l.[56j. 

By way of contrast, RCC matrix element ratios ex- 
hibit about 4% differences from the geometric ratios. One 
would expect the RCC matrix element ratios to be much 
closer to the MBPT-SD ratio given the close formal sim- 
ilarities between the RCC and MBPT-SD approaches. 



The RCC matrix element ratios listed in Table IIIII also 
show significant differences from the geometric ratio for 
the 4s -^ 3dj transitions. The DFCP and MBPT-SD 
ratios lie within 1% of the geometric ratios. It should be 
noted that a similar situation exists for the 5s — 4^5/2:3/2 
line strength ratios of Sr+ with RCC calculations exhibit- 
ing much larger differences due to non-geometric effects 
than other calculations ;57] . The feature common to the 
DFCP and MBPT-SD methods is that they use large B- 
spline basis sets and calculated quantities are expected 
to be independent of basis set effects. One possible cause 
for the different RCC matrix element ratios lies in the 
gaussian basis set used to represent virtual excitations in 
the RCC calculation. This point will be addressed later 
where polarizabilities are discussed. 



TABLE IV: Lifetime of the 8^3/2 and 8^5/2 levels of Ca+ (in 
sec). The 8^3/2 : 8d5/2 lifetime ratio is also given. 



Source 



T"3d 



DFCP 

MBPT-SD [29] 
RCC [26] 
MCHF [34] 
Experiment [23] 
Experiment [58] 
Experiment [59] 
Experiment [60] 



3/2 



T-idr, 



L143(l)(s) 

L196(l)(s) 

1.185(7) 

1.160 
1.176(11) 

1.17(5) 
1.064(17) 
1.111(46) 



Z2_ 



Ratio 



1.114(1) 
1.165(11) 
1.110(9) 

1.140 
1.168(7) 
1.09(5) 



1.0260 
1.0266 
1.0675 
1.0175 
1.007(15) 
1.073(90) 



0.994(38) 1.118(80) 



TABLE V: Lifetimes (in nsec.) of the 4pj. and 4p3 states. 
The 4pi : 4p3 lifetime ratio is also given. The quantity R 
gives fraction of the total decay rate arising from the indicated 
transition. 



Level DFCP MBPT-SD RCC 

[29] [26] 



4pi (ns) 
4p3 (ns) 



7?(4pi -4si; 
i?(4pl -3dl] 
i?(4p3 - 4s 1) 
_R(4p3 -3dl] 
Riipl ~3dl] 
Ratio 



6.94(1) 
6.75(1) 



0.9324 
0.0676 
0.9313 
0.0069 
0.0617 
1.0281 



6.88(6) 



6.931 



1.69(6) 6.881 



0.9340 
0.00667 
0.0593 
1.0284 



0.9374(74) 

0.0626(5) 

0.9350(62) 

0.00666(4) 

0.0583(4) 

1.0073 



Expt. 

7.098(20) [66] 
7.07(7) [67] 
6.926(19) [66] 
6.87(6) [671 
6.72(2) [68| 
6.61(30) [69| 



0.9347(3) [70] 
0.00661(4) [70] 
0.0587(2) [70] 
1.025(3) [66] 
1.029(14) [67] 



B. Lifetimes 



lifetimes closer to the MBPT-SD and DFCP hfetimes. 

Measurements of the branching ratios of the 4p3/2 state 
yield a picture where the MBPT-SD calculations largely 
agree with experiment while the DFCP tends to overes- 
timate the contributions of the decays to the 3dj levels. 
Another area of partial agreement between theory and 
experimental occurs for the 4px/2 ■ 4^3/2 lifetime ratio. 
The DFCP, MBPT-SD and experimental ratios range 
from 1.025 to 1.030, with the RCC calculation again pro- 
viding an outlier at 1.0073. 



The two most important lifetimes for the Ca^ clock 
[6ll - l65l ] are the lifetimes of the 3dj and 4pj levels. 
The 3dj states decay to the ground state in an electric 
quadrupole transition with lifetime of about 1.1 sec [23| . 
The 4pj states experience electric dipole transitions to 
both the 3d J and 4s states. Table ITVl gives the lifetimes 
of the 3d J states while Table |V] gives the lifetimes of the 
4pj states. All DFCP lifetimes were computed using ex- 
perimental energy differences. 

The most recent experiment for the 3dj lifetimes give 
a ratio of 1.007 ±0.015 sec for the 3^3/2 and 3^5/2 states. 
This suggests that the 4s -^ 3dj matrix element ra- 
tios should be close to the values expected from angular 
momentum coupling considerations. Older experiments 
[53, |60] give ratios further from unity, but in these cases 
the uncertainties are much larger. 

The lifetimes of the Apj states depend on two transi- 
tions, these are the 4s-4pj and 3dji-Apj transitions, with 
the 4s-4pj transition being the most important. The life- 
times and branching ratios for the Apj states are given in 
Table FVl It not possible to reconcile the theoretical and 
experimental lifetimes at the 1% level. The two most re- 
cent experiments 66, 67] gave lifetimes that are 2% larger 
than the DFCP lifetimes and 3% larger than the MBPT- 
SD lifetimes. Older Hanle effect experiments |6a,i62i] gave 



IV. POLARIZABILITIES 
A. Static Polarizabilities 

The static dipole and quadrupole polarizabilities are 
calculated by the usual sum-rule 



a 



W 



E 



/, 



2 



(11) 



pOT 



where the /^^ are the absorption oscillator strengths and 
Egi is the excitation energy of the transition. Static dipole 
polarizabilities for the 4s, 4pj and 3dj states are listed 
in Table I VI I All polarizabilities were computed using ex- 
perimental energy differences. 

The most important polarizability is that of the 4s 
ground state and there is only a 1% variation between 
the DFCP, MBPT-SD and CICP static dipole polariz- 
abilities. The DFCP polarizability is smaller than the 
MBPT-SD polarizability because the DFCP 4s - Apj 
matrix elements are smaller. The RCC calculation of 
the dipole polarizability is the clear outlier at 73.0 a.u. 
[26t . The good agreement between the DFCP, CICP and 



TABLE VI: Dipole and quadrupole polarizabilities (in a.u.) for low-lying states of the Ca^ ion. 
polarizabilities are not given for states with £ > 0. The RCC-STO results are those from Ref. 
orbital basis to represent virtual excitations. 



Non-relativistic quadrupole 
i26j that used a Slater type 



State DFCP 



„(0) 



v(*) 



Others 



DFCP 



Others 



DFCP 



Others 



4si/2 75.28 

4pi/2 -2.774 
4p3/2 -0.931 



34/2 32.99 



3d5/2 32.81 



76.1(5) MBPT-SD [29] 

75.49 CICP [24] 

73.0(1.5) RCC [26] 

75.3(4) /-sums [72] 

74.3 RCC-STO [26] 

-0.75(70) MBPT-SD [29] 

-2.032 CICP [2^ 
1.02(64) MBPT-SD \29} 

-2.032 CICP \2£ 

32.0(3) MBPT-SD [29] 

32.73 CICP [24] 

28.5(1.0) RCC [26] 

31.6 RCC-STO [26] 

31.8(3) MBPT-SD [29] 

32.73 CICP [2J 

29.5(1.0) RCC [26] 

32.5 RCC-STO [261 



10.12 



-17.1 



-25.16 



10.31(28) MBPT-SD [29] 

10.47 CICP [24] 

-17.43(23) MBPT-SD [29] 

-17.64 CICP [24] 

-15.87 RCC [26] 

-17.7 RCC-STO [26] 

-24.51(29) MBPT-SD [29| 

-25.20 CICP \2£ 

-22.49(5) RCC [26] 

-25.5 RCC-STO [26] 



882.43 



7.466 [4] 

-3.571 [4] 

4928 



-3304 



871(4) MBPT-SD [29] 

875.1 CICP [24] 

712.9(24) RCC [71] 

906(5) RCC [45] 



-3706(75) RCC [45] 



MBPT-SD polarizabilities does not necessarily imply a 
1% reliability in these polarizabilities since the calcu- 
lations give lifetimes for the Apj states that are 2-3% 
smaller than experiment. 



The variation between the DFCP, MBPT-SD and 
CICP estimates of the 3dj state polarizabilities do not 
exceed 1.0 a.u. The difference in the polarizabilities for 
the two members of the spin-orbit doublet is only 0.2 a.u. 



The polarizabilities of the Apj states are close to zero 
with the polarizability of the 4p3/2 state being about 1.8 
a.u. larger than the polarizability of the 4pi/2 state. 
The polarizability is small because the downward transi- 
tions to the 4si /2 and 3dj states have negative oscillator 
strengths which result in cancellations in the oscillator 
strength sum. This is evident in Tables I VIII and IVIIII 
which show the breakdown of the different contributions 
to the polarizabilities from the oscillator strength sum 
rule. 



TABLE VII: The contributions of individual transitions to 
the polarizabilities of the 4si/2 and 4pi/2 states at the magic 
wavelengths. The numbers in brackets are uncertainties in 
the last digits of the energy or wavelength calculated by in- 
troducing 2% uncertainties into the most important matrix 
elements. 



cj (a.u.) 0.0659561(11247) 

A (nm) oo 690.817(11.984) 



0.1152981(4) 0.1238091(303) 
395.1807(14) 368.0149(901) 







4Si/2 






4pi/2 


24.0704 


35.9364 


-2665.2940 


-147.2228 


5pi/2 


0.0097 


0.0102 


0.0117 


0.0121 


4p3/2 


47.7532 


70.6856 


5558.6017 


-333.5265 


5P3/2 


0.0145 


0.0153 


0.0175 


0.0181 


Remainder 


0.1672 


0.1710 


0.1794 


0.1815 


Core 


3.2600 


3.2664 


3.2793 


3.2823 


Total 


75.2751 


110.0849 

4pi/2 


2896.7954 


-477.2554 


4Si/2 


-24.0704 


-35.9364 


2665.2940 


147.2228 


5Si/2 


11.7449 


16.4949 


97.8655 


-798.9861 


3(^3/2 


-39.6152 


69.1092 


10.4051 


8.7196 


4^3/2 


40.8730 


51.6866 


113.3267 


155.5319 


Remainder 


5.0332 


5.4542 


6.6245 


6.9740 


Core 


3.2600 


3.2664 


3.2793 


3.2823 


Total 


-2.7742 


110.0850 


2896.7954 


-477.2554 



The comparisons of the polarizabilities sug gest that 
the basis set used in the RCC calculations [261] could be 
improved. The recommended results for the RCC calcu- 
lation are those computed with the gaussian basis. How- 
ever, RCC calculations performed using a Slater type 
orbital basis [2^ give polarizabilities that are in much 
better agreement with the MBPT-SD and DFCP polar- 
izabilities. 



B. Dynamic polarizabilities and magic w^avelengths 

The dynamic dipole polarizability of a state at photon 
energy to is defined 



ai 



H-E 



^(1) 

I gt 



£^ ■ —Up' 
9* 



(12) 



TABLE VIII: The contributions of individual transitions to the polarizabihties of the 4si/2 and 4p3/2 states at the magic 
wavelengths. These results assume non-polarized light. The numbers in brackets are uncertainties in the last digits calculated 
by assuming certain matrix elements have ±2% uncertainties. 



iu (a.u.) 





0.0663204(11651) 


0.1149923(4) 


0.1232650(511) 


0.0677517(11210) 


0.1151251(3) 


A (nm) 


CXD 


687.022(12.285) 


396.2315(13) 


369.6393(1534) 


672.508(11.3150) 


395.7748(10) 


4pi/2 


24.0704 


36.1337 


-6530.5659 


4Si/2 

-157.0218 


36.9414 


-4009.0830 


5pi/2 


0.0097 


0.0103 


0.0117 


0.0121 


0.0103 


0.0117 


4p3/2 


47.7532 


71.0638 


3449.6093 


-358.6372 


72.6100 


4129.0858 


5p3/2 


0.0145 


0.0154 


0.0175 


0.0181 


0.0154 


0.0175 


Remainder 


0.1672 


0.1710 


0.1794 


0.1813 


0.1712 


0.1794 


Core 


3.2600 


3.2664 


3.2791 


3.2820 


3.2667 


3.2792 


Total 


75.2751 


110.6606 


-3077.3881 


-512.1655 

4p3/2 

rUj = 1/2 


113.0150 


123.4906 




Average 


m, = 1/2 


rrij = 1/2 


ruj = 3/2 


rrij = 3/2 


4Si/2 


-11.9383 


-71.0636 


-3449.6902 


358.6371 


0.0000 


0.0000 


5Si/2 


6.0501 


34.3769 


219.9519 


-1069.3049 


0.0000 


0.0000 


3^3/2 


- 5.4283 


1.4608 


0.2153 


0.1808 


11.6778 


1.9317 


4^3/2 


5.8429 


1.0626 


2.3512 


3.2245 


9.6796 


21.2507 


3^5/2 


-31.6965 


77.2921 


11.5779 


9.7303 


45.8671 


7.6960 


44/2 


33.7190 


57.2196 


126.3638 


173.0763 


38.6057 


84.5964 


Remainder 


4.3193 


7.0456 


8.5623 


9.0081 


3.9179 


4.7365 


Core 


3.2600 


3.2664 


3.2791 


3.2820 


3.2667 


3.2792 


Total 


- 4.1279 


110.6606 


-3077.3881 


-512.1655 


113.0150 


123.4906 
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FIG. 1: (color online) Dynamic polarizabilities of the 4si/2 
and 4pi/2 states of the Ca"*" ions. Magic wavelengths are 
identified by arrows. 



The dipole polarizability has a tensor component for 
states with states with J > 1/2. This can be written 



a^ (uj) — 6 
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TABLE IX: Pseudo-spectral oscillator strength distribution 
used in the computation of the dynamic polarizability of the 
Ca^"*" core. Energies are given in a.u.. 



133.689002 
14.645933 
11.675258 
1.9047772 
1.1104171 



2.0 
2.0 
6.0 
2.0 
6.0 



The polarizability for a state with non-zero angular mo- 
mentum J depends on the magnetic projection Mg: 



aiM„ 
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The dynamic polarizabilities includes contributions 
from the core which is represented by a pseudo-oscillator 
strength distribution [3ll[73,[7J] which is tabulated in Ta- 
ble llXI The distribution is derived from the single particle 
energies of a Hartree-Fock core. Each separate (n, t) level 
is identified with one transition with a pseudo-oscillator 
strength equal to the number of electrons in the shell. 
The excitation energy is set by adding a constant to the 
Koopmans energies and adjusting the constant until the 
core polarizability from the oscillator strength sum rule 
is equal to the known core polarizability of 3.26 a.u. The 
core polarizabilities of any two states effectively cancel 
each other when the polarizability differences are com- 
puted. 

The dynamic polarizabilities for the 4si/2 and 4pi/2 
states of Ca"*" are shown in Figure [T] The first magic 
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FIG. 2; (color online) Dynamic polarizabilities of the 4si/2 
and 4p3/2 states of Ca^ . Magic wavelengths are identified by 
arrows. 
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FIG. 3; (color online) Dynamic polarizabilities of the 4si/2 
and 3(^5/2 states of Ca^ . Magic wavelengths are identified by 



wavelength occurs a.tu! — 0.0659561 a.u. after the photon 
wavelength exceeds the energy for the 4^1/2-3^3/2 tran- 
sition. Magic wavelengths are identified at A = 690.817, 
395.181 and 368.015 nm. The 395.181 nm magic wave- 
lengths occur when the photon is very close to the excita- 
tion energies of the 4s-4pj states. The 368.015 nm magic 
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FIG. 4: (color online) Dynamic polarizabilities of the 4si/2 
and 3d5/2 states of Ca"^. Magic wavelengths is identified by 
circles and arrows. 



wavelength occurs near the energy for the 4pi/2 — 5si/2 
transition. The dominant contributions to polarizabili- 
ties at the magic wavelengths are listed in Table Ivnl The 
4s polarizability is dominated by the 4s]^/2-4pj transi- 
tions with the next largest contribution coming from the 
core. However, the 4pi/2 polarizability has significant 
contributions from the transitions to the 4s, 5s and 3^3/2 
states. A magic wavelength experiment would give infor- 
mation about the 4pi/2 state, but would not give detailed 
information about any individual matrix element. An 
experiment that measured all three magic wavelengths 
could conceivably be able to extract information about 
individual line strengths, however it should be noted that 
two of the transitions are in the ultraviolet. 

The dynamic polarizabilities of the 4si/2 and 4p3/2 
states of Ca"*" are shown in Figure O These figures as- 
sume non-polarized light. Figure [2] only has two magic 
wavelengths below w = 0.125 a.u. Transitions to the 
nsi/2 states make no contribution to the 4p3/2 state po- 
larizability. This is evident from Table I Villi which de- 
tails the breakdown of different transitions to the polar- 
izability. The magic wavelength at 395.775 nm for the 
4p3/2,m=3/2 niagnetic sub-level can give an estimate of 
the contribution to the np-j^i^ polarizability arising from 
excitations to the ndj levels. 

The 4si/2 and 3d5/2 polarizabilities are shown in Fig- 
ures [3] and m The 3(i5/2_,„ polarizabilities are shown for 
all magnetic sub-levels and also for the average polariz- 
ability. Magic wavelengths occur when the photon energy 
gets close to the excitation energies for the 3^5/2 -^ 4p,/ 
transitions and the 4si/2 — >■ 4pj transitions. Figure [3] 
shows the 4si/2 and 3^5/2 polarizabilities at photon en- 
ergies between 0.02 and 0.07 a.u. Precise values of the 
magic wavelengths and the breakdown of the polarizabil- 
ity into different components can be found in Table [X] 

Two of the magnetic sub-levels have magic wave- 
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lengths at infrared frequencies, namely A = 1338.474 
and 1074.336 nm. The contributions to the in 3^5/2 po- 
larizability are dominated by the 3^5/2 — ^ '^Pi/2 transi- 
tion which constitutes about 88% of the polarizability. 
The measurement of these magic wavelengths provides a 



method to determine the // 



4si/2- 



to /sds. 



P3/2 



oscil- 



lator strength ratios. Suppose all the remaining compo- 
nents of the 3(i5/2 polarizability can only be estimated to 
an accuracy of 10%. The overall net uncertainty in the 
remaining terms would be less than 1.5%. 

There are also an additional magic wavelengths that 
can potentially be measured. The 4s dynamic polar- 
izability goes through zero as the wavelength passes 
through energies needed to excite the 4s — > 4pi/2 and 
4s — >■ 4j>3/2 transitions. Figure 2] shows the polarizabili- 
ties for the 4s and 3^5/2 at energies near the 4s — ?> 4pj 
excitation energies. The 3^5/2 polarizabilities are typi- 
cally small in magnitude in this wavelength range. The 
magic wavelength arises more from the the cancellation 
of the 4pi/2 and 4p3/2 contributions to the 4s dynamic 
polarizability than from the cancellation between the 4s 
and 3^5/2 dynamic polarizabilities. Measurement of the 
magic wavelength here is in some respects in analogous 
to a measurement of the longest tune-out wavelength 
for neutral potassium [56]. Zero field shift wavelengths 
measured in the spin-orbit energy gap of the resonant 
transition are strongly dominated by the large and op- 
posite polarizability contributions of the two members of 
the spin-orbit doublet [5a, uM- This makes it possible 
to accurately determine the oscillator strength ratio, i.e. 
f4:s^Api,2 ■/4s->4p3/2 1 of the two transitions comprising the 
spin-orbit doublet. 
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0.04 0.05 
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FIG. 5: (color online) Dynamic polarizabilities of the 4si/2 
and 3(^3/2 states of Ca"*". Magic wavelengths are identified by 
arrows. 

Table IXII identifies the magic wavelengths associated 
with the 4s — > 3^3/2 energy interval. The situation here 
is similar to the situation for the 4s — >■ 3^5/2 magic wave- 



lengths. However, there are three magic wavelengths in 
the infrared region of the spectrum. This transition has 
an additional magic wavelength since the 3(i3/2,m=i/2 
state, unlike the 3^5/2, m=i/2 state, also undergoes un- 
dergoes a transition to the 4:pi/2 state. The polarizabil- 
ity difference in the 0.02 to 0.07 a.u. energy range is 
plotted in Figure [5l The 3^3/2 polarizability is domi- 
nated by the 3^3/2 -^ 4pj transition and a magic wave- 
length measurement can be used to make an estimate 
of the 3(^3/2 -^ ipj line strength relative to the 4s dy- 
namic polarizability. The 3^3/2, m=i/2 polarizability at 
850.335 nm has large contributions from the 4p]^/2 and 
4p3/2 states since it lies between the excitation ener- 
gies of these of states. Measurement of the 850.335 
nm and 1308.590 nm wavelengths together would give 



estimates of the 3d 



fsd^ 



■fsds 



3/2 ~^ '^Pi/2 liii6 strengths and the 
ratio. A measurement of the 



^3/2^4^1/2 •./3d3/2^4p3/2 

magic wavelengths in the vicinity 395 nm provides would 
permit a determination of the fis^4pi/2-f4.s-y4p3/2 ratio. 



C. Uncertainties 

An uncertainty analysis has been done for all the magic 
wavelengths presented in the preceding sections. This 
analysis was aimed at making an initial estimate of how 
uncertainties in the matrix elements of the most impor- 
tant transitions would translate to a shift in the magic 
wavelengths. The primary purpose of the uncertainty 
analysis is to define reasonable limits to help guide an 
experimental search for the magic wavelengths identified 
in this paper. 

In the case of the 4s -> Apj polarizability differences, 
the 4s -^ Apj, Apj -^ 5s, Apj -^ 3dj and Apj -^ Adj 
matrix elements were all changed by 2% and the magic 
wavelengths recomputed. The matrix elements involving 
the different spin-orbit states of the same multiplet were 
all given the same scaling. A variation of ±2% was cho- 
sen by reference to the difference of the DFCP matrix 
elements with the experimental or the MBPT-SD ma- 
trix elements. The estimate of a 2% uncertainty in the 
4s -^ Apj matrix element can be regarded as a conserva- 
tive estimate. 

The 4s -^ 3dj polarizability difference is predomi- 
nantly determined by the 4s — s- 4p,/ and 3dj — >■ Apj 
matrix elements. So variations of ±2% in these two tran- 
sitions were used in determining the uncertainties in the 
magic wavelengths. 

There are a number of magic wavelengths which are 
relatively insensitive to changes in the matrix elements 
of a multiplet. One of these wavelengths is the 850 nm 
wavelength for the 4s — 3^3/2 interval and the others are 
the magic wavelengths near 395 nm. These wavelengths 
arise due to cancellations in the polarizabilities due to 
two transitions of a spin-orbit doublet. In the case of the 
850 nm magic wavelength, the relevant transitions are 



the 3d' 



3/2 - 



►4pj transitions. 



The sensitivity of the magic wavelengths near 395 nm 



to changes in the transition matrix elements depends on 
the overall size of the polarizabilities of the Apj and 3dj 
levels. When these are large due to transitions other 
than the 4s — > Apj transition, then the 395 nm magic 
wavelength shows higher sensitivity to the changes in the 
matrix elements. However, the net change in the magic 
wavelengths for 2% changes in the matrix elements is 
about 0.001 nm for the 4s — )► 4p,/ interval. The sensi- 
tivity to 2% matrix element changes for the 4s — > 3dj 
intervals is about 0.0001 nm due to the small polarizabil- 
ities of the 3dj states near 395 nm. The 850 nm magic 
wavelength is also relatively insensitive to changes in the 
overall size of the matrix elements, with the 2% matrix 
element change leading to a change of only 0.0001 nm 
in the magic wavelengths. The low sensitivity of magic 
wavelengths to the overall size of the matrix elements in 
these cases means that these the magic wavelengths can 
be used to give precise estimates of the matrix element 
ratios of the two transitions in the spin-orbit doublet. 

The 1338, 1309, 1074, 887 nm magic wavelengths show 
much greater sensitivity to 2% changes in the matrix ele- 
ments. The changes in the magic wavelengths range from 
3 to 80 nm. The sensitivity of the magic wavelengths to 
these matrix elements is driven by the rate of change 
of the 4s and 3dj polarizabilities with energy. A large 
change in the photon energy is needed to compensate for 
a small change in the polarizability when dai/dio is small. 
The sensitivity of the magic wavelength to small changes 
in the matrix elements decreases as the photon energy 
gets closer to the 3dj — J> 4pj/ excitation thresholds. The 
high sensitivity of the magic wavelengths with respect to 
changes in the matrix elements means it is only necessary 
to measure the magic wavelength to a precision of 0.10 
nm to impose reasonably tight constraints on the ratios 
of the 4s— >-4pj and 3(ij— >4pj' matrix element rations. 



sonably accurate estimates of the Sdj-Apj' line strengths 
relative to the As-ipj line strengths. This could lead 
to improved estimates of the blackbody radiation shift 
for the Ca"^ clock transition. There is one impediment. 
At the moment there is a 3% spread between theoreti- 
cal and experimental lifetimes for the 4p,/' states. This 
variation, which does not exist for the same transition in 
potassium [56, 76], needs to resolved so the uncertainty 
in the 4s — Apj line strengths can be reduced to 1% or 
better. 

There are two other relatively clean measurements of 
atomic structure parameters that could be made. Mea- 
surement of the magic wavelength near 395 nm could 
be used to determine a value of the oscillator strength 
f4s-^4pi,2-f4s^4p-i,2 ratio. This could help resolve the in- 
compatible predictions of this ratio by DFCP/MBPT- 
SD and RCC calculations. Comparisons of polarizabili- 
ties do suggest that the gaussian basis set used for the 
RCC calculations could be improved. Further, mea- 
surements of the two longest magic wavelengths for the 
3^3/2. m=i/2 "^ 4si/2 transition could give a good esti- 
mate of the /3d3/2^4pi/2:/3d3/2^4p3/2 r^tio. 

The utility of measuring magic wavelengths for se- 
lected Ca+ transitions can of course be extended to other 
alkaline-earth ions, with Sr^" and Ba+ being obvious pos- 
sibilities. A single ion optical frequency standard at the 
10~^^ level of precision has recently been reported for 
the 5s-4(i5/2 transition of the Sr"*" ion 77[. It is likely 
that the determination of the magic wavelengths for this 
transition could be used to improve the precision of esti- 
mates of the blackbody radiation shift for this transition 

MM. 
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TABLE X: The contributions of individual transitions to the polarizabilities of the 4si/2 and 80(5/2 states at the magic 
wavelengths. These results assume non-polarized light. The numbers in brackets are uncertainties in the last digits calculated 
by assuming certain matrix elements have ±2% uncertainties as described in the text. 



u) (a.u.) 





0.0340414(22387) 


0.0424109(10654) 


0.1151182(1) 


0.1151184(1) 


0.1151186(1) 


A(nm) 


CXD 


1338.474(82.593) 


1074.336(26.352) 


395.7982(1) 


395.7978(1) 


395.7968(1) 


4pi/2 


24.0704 


26.3917 


27.8762 


4Si/2 


-4090.5249 


-4088.7574 


-4085.2247 


5pi/2 


0.0097 


0.0098 


0.0099 




0.0117 


0.0117 


0.0117 


4p3/2 


47.7532 


52.2705 


55.1513 




4087.5752 


4088.4488 


4090.2003 


5p3/2 


0.0145 


0.0147 


0.0148 




0.0175 


0.0175 


0.0175 


Remainder 


0.1672 


0.1682 


0.1688 




0.1794 


0.1794 


0.1793 


Core 


3.2600 


3.2618 


3.2627 




3.2792 


3.2792 


3.2792 


Total 


75.2751 


82.1167 


86.4837 


3^5/2 


0.5371 


3.1792 


8.4633 




Average 


mj = 1/2 


ruj = 3/2 


m, = 1/2 


rrij = 3/2 


rrij = 5/2 


4p3/2 


29.5834 


71.3309 


76.6749 




-11.5457 


-7.6971 


0.0000 


5p3/2 


0.0113 


0.0165 


0.0119 




0.0227 


0.0151 


0.0000 


4/5/2 


0.0607 


0.0109 


0.0988 




0.0136 


0.1223 


0.3398 


5/5/2 


0.0196 


0.0035 


0.0318 




0.0041 


0.0367 


0.1018 


4/7/2 


2.5573 


3.2582 


2.7444 




4.0780 


3.3983 


2.0391 


5/7/2 


0.8270 


1.0479 


0.8799 




1.2223 


1.0186 


0.6112 


Remainder 


2.5979 


3.1870 


2.7803 




3.4628 


3.0060 


2.0922 


Core 


3.2600 


3.2618 


3.2627 




3.2792 


3.2792 


3.2792 


Total 


38.5915 


82.1167 


86.4837 




0.5371 


3.1792 


8.4633 



TABLE XL The contributions of individual transitions to the polarizabilities of the 4si/2 and 3d3/2 states at the magic 
wavelengths. These results assume non-polarized light. The numbers in brackets are uncertainties in the last digits calculated 
by assuming certain matrix elements have ±2% uncertainties as described in the text. 



U) 





0.0348188(20007) 


0.0513460(1855) 


0.0535831(1) 


0.1151182(1) 


0.1151185(1) 


A 


00 


1308.590(71.108) 


887.382(3.196) 


850.335(2) 


395.7981(1) 


395.7970(1) 


4pi/2 


24.0704 


26.5098 


4Sl/2 

30.0922 


30.7777 


-4090.1007 


-4086.1347 


5pi/2 


0.0097 


0.0098 


0.0100 


0.0101 


0.0117 


0.0117 


4p3/2 


47.7532 


52.4999 


59.4402 


60.7642 


4087.7840 


4089.7488 


5p3/2 


0.0145 


0.0147 


0.0150 


0.0151 


0.0175 


0.0175 


Remainder 


0.1672 


0.1683 


0.1695 


0.1697 


0.1794 


0.1794 


core 


3.2600 


3.2619 


3.2639 


3.2643 


3.2792 


3.2792 


Total 


75.2751 


82.4644 


92.9908 

3d3/2 

rrij = 3/2 


95.0011 


1.1711 


7.1019 




Average 


rrij = 1/2 


rrij = 1/2 


m, = 1/2 


rrij = 3/2 


4pi/2 


9.9038 


70.5419 





-1034.8996 


-10.4461 





5pi/2 


0.0033 


0.0134 





0.0139 


0.0184 





4p3/2 


5.4284 


1.3416 


84.7097 


1119.1728 


-0.2147 


-1.9320 


5p3/2 


0.0021 


0.0003 


0.0029 


0.0003 


0.0004 


0.0038 


4/5/2 


2.3339 


3.1745 


2.1676 


3.2644 


3.9674 


2.6449 


5/5/2 


0.7556 


1.0218 


0.6928 


1.0422 


1.1908 


0.7939 


Remainder 


2.3601 


3.1091 


2.1539 


3.1428 


3.3757 


2.3121 


core 


3.2600 


3.2619 


3.2639 


3.2643 


3.2792 


3.2792 


Total 


24.0472 


82.4644 


92.9908 


95.0011 


1.1711 


7.1019 



